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ABSTRACT. Study of the polarization of supernovae has suggested that the core collapse process 
may be intrinsically strongly asymmetric. There is a tentative trend for supernova with smaller en- 
velopes showing more polarization, with SN Ic having the smallest envelopes and showing the largest 
polarization. The recent discovery of the unusual supernova SN 1998bw and its apparent correlation 
with the gamma-ray burst GRB 980425 has raised new issues concerning both the 7-ray bursts and 
supernovae. SN 1998bw resembled a SN Ic, but was unusually bright at maximum light in the optical 
and radio, and its expansion velocities were large. This makes SN 1998bw a possible candidate for a 
"hypernova" with explosion energies exceeding 10 52 erg. We show that the light curve of SN1998bw 
can be understood as the result viewing an aspherical explosion roughly along the symmetry axis of 
an exploding, non-degenerate C/O core of a massive star with a kinetic energy of 2 X 10 51 erg, a total 
ejecta mass of 2Mq, and a 56 7Vi mass of 0.2Mq. In this model, the high expansion velocities are a 
direct consequence of the aspherical explosion which, in turn, produces oblate iso-density contours and 
that accounts for the polarization. It is not yet clear how either the hypernovae or these asymmetric 
models can produce 7-ray bursts. 

1. Introduction 

Due to its correlation in time and location, the 7-ray burst GRB 980425 has a high 
probability of being associated with SN 1998bw (Galama et al. 1998). This connection 
is supported by the association of a radio source with SN 1998bw that requires rapid 
expansion (Waxman & Loeb 1998) and probably relativistic expansion (Kulkarni et al. 
1998). Recent corrections to the BeppoSAX positions and record of time variability 
show that the supernova falls in the error box (Piro, et al. 1998) of a time-variable 
X-ray source (Pian, et al. 1998). From optical spectra, SN 1998bw was classified as a 
SN Ic by Patat and Piemonte (1998). What sets SN 1998bw apart from other Type 
Ic (SN Ic) are higher expansion velocities as indicated by the Si II and Ca H and K 
lines (« 30 to 50% higher at maximum light than SN1994I and SN1983V; Clochiatti & 
Wheeler 1997), the red colors at maximum light, and the large intrinsic brightness. A 
peak luminosity of 1.3 ± 0.6 x 10 43 erg s _1 can be inferred from the redshift of the host 
galaxy (z = 0.0085, Tinney et al. 1998) and the reddening (A v = 0.2 m , Schlegel et al. 
1998), if we assume H a = 67 km/s/Mpc. The uncertainties are rather large as the host 
galaxy is not yet fully in the Hubble flow, so the peculiar velocity may be of the order 
of 300 to 400 km s _1 , H a is known only to an accuracy of ss 10%, and Ay may vary by 
0.1 m . The radio source associated with SN 1998bw is the brightest ever observed for 
a supernova (Kulkarni, et al. 1998). 



The properties of SN 1998bw suggest that it was a "hypernova" event (Paczyhski 
1997). Based on their light curve calculations, Iwamoto et al. (1998) and Woosley, East- 
man & Schmidt (1998) d erived explosion energies of 20-50 foe and 22 foe (1/oe = 
10 51 erg) ejecta masses of 12-15 M Q and 6 M , and 56 Ni masses of 0.6-0.8 and 0.5 M Q , 
respectively. 

Guided by the deduced properties of more traditional core collapse supernovae and 
SN Ic in particular, we have investigated whether asphericity can provide an alternative 
explanation for SN 1998bw that allows its optical properties, at least, to remain in the 
range of "normal" SN Ic. 

2. Polarization of Supernovae 

Both spectral analyses and light curve calculations support the picture that SN II, SN lb 
and SN Ic form a sequence involving core collapse with successively smaller H and He 
envelopes (Clochatti & Wheeler 1997). The analysis of spectra and light curves gives, 
however, essentially no insight into the geometry of the expanding envelope. Polarization, 
on the other hand, provides a unique tool to explore asymmetries. For the last several 
years we have engaged in a program to obtain spectropolarimetry of as many supernovae 
as possible from McDonald Observatory (Wang et al. 1996). Linear polarization of « 1% 
seems to be typical for SN II (Wang, Wheeler & Hoflich 1998). There is a trend, however, 
for the observed polarization to increase in core-collapse supernovae with decreasing 
envelope mass, e.g. from SN II to SN Ic (Wang et al. 1998). SN 1987A, with a hydrogen 
envelope of about 1OM , had a polarization of w 0.5%; for SN 1993 J, with a hydrogen 
envelope mass of only a few tenths of M Q , the observed linear polarization was as high 
as « 1.5%; for the SN Ic 1997X, with no hydrogen and little or no helium envelope, 
the polarization was even higher, perhaps several percent (Wang et al. 1998; Wang & 
Wheeler 1998). The suggestion is that the closer one looks to the collapsing core itself, 
the larger the polarization and the larger the asymmetry. This trend, while tentative, 
clearly points toward the interpretation that the explosion itself is strongly asymmetric. 

From theoretical calculations of scattering dominated atmospheres, this size of po- 
larization, <; 1 %, requires axis ratios of the order of 2 or 3 to 1, requiring the ejecta 
to be highly aspherical. The luminosity L(Q) will vary by a factor of « 2 as the line of 
sight varies from the equator to the pole (Hoflich 1991, Hoflich et al. 1995). 

Given the ubiquitous presence of polarization in core collapse supernovae and es- 
pecially SN Ic, inclusion of asphericity effects in SN Ib/c may prove to be critical to 
their understanding (Wang et al. 1998). Polarization was been observed in SN1998bw 
(Kay et al. 1998) at the level of about 0.5%. Asymmetries thus cannot be neglected in 
a complete model for this event. An obvious question is whether the asymmetries are 
directly connected to the fact that it produced a 7-ray burst. 

3. Type Ib/c Supernovae and Gamma- Ray Bursts 

Wang & Wheeler (1998) examined the question of whether there is a general connection 
between supernovae and 7-ray bursts. They concluded that SN la could be excluded from 
any such correlation at the 4cr level, that the sample of sufficiently well-studied SN II 



was too small to reach meaningful conclusions, but that a correlation with SN Ib/c could 
not be ruled out. Table I gives the list of candidate SN Ib/c - 7-ray burst correlations 
suggested by Wang & Wheeler. 



Table 1 - Supernovae of Type Ib/c and GRB Associations 



SN 


Type 


Dates 


RA (2000.0) 


DEC (2000.0) 


SN 1994at 


Ib/c 


-961009 


17.1 


-1.0 


GRB 960925 






29.37 


-13.9 


SN 1996N 


lb 


-960310 


54.7 


-26.3 


GRB 960221 






47.8 


-31.24 


SN 1997B 


Ic 


-970104 


88.3 


-17.9 


GRB 971218 






97.75 


-21.73 


SN 1997dq 


lb 


-971105 


175.2 


+11.5 


GRB 971013 






167.0 


+2.7 


SN 1997ef 




-971205 


119.3 


+49.6 


GRB 971152 






84.6 


+41.7 


SN 1997ei 


Ic 


-971223 


178.5 


+58.5 


GRB 971120 






155.7 


+76.4 


SN 1998bw 


7 


980424-980427 


19 35 03.3 


-52 50 44.8 


GRB 980425 




980425.909 


19 35 21 


-52 52 19 



Kippen et al. (1998) found no evidence for a connection between supernovae and 
7-ray bursts and Graziani, Lamb, & Marion (1998) concluded that not all SN Ib/c 
could be associated with 7-ray bursts, although a fraction might be. With the small 
sample of SN Ib/c and partial sky coverage, it will be difficult to resolve this issue 
with statistics alone. Woosely, Eastman, & Schmidt (1998) suggested that SN 1998cy 
might be associated with GRB 9970514. Unlike the suggestions in Table I, SN 1998cy 
shows both narrow and broad features of hydrogen. This event was even brighter than 
SN 1998bw (Schmidt, 1998). If it is correlated with a 7-ray burst it represents an- 
other type of exploding object. When it was first observed, SN 1997ef was noted to 
be an object of unprecedented spectral features, although it bore some resemblance to 
SN Ic. The similarity of SN 1997ef to SN 1998bw was pointed out by Garnavich (1997). 
Wang & Wheeler found that SN 1997ef was just beyond the 2cr BATSE error box for 
GRB 971115. A model for the light curve of SN 1997ef has been presented by Iwamoto 
et al. (1999) based on a carbon/oxygen core exploding with a "normal" kinetic energy. 
At this meeting, Nomoto reported a revised calculation based on a "hypernova" scenario 
with in excess of 10 foe of kinetic energy that better matched the line profiles. The ques- 
tion of whether this event was polarized and hence asymmetric also arises. Garnavich, 
Jha, & Kirshncr (1998) report that SN 1998ey was "identical" to SN 19978ef two weeks 
after discovery and similar to SN 1998bw. They report no correlation with a BATSE 
trigger, but again, BATSE only surveys a portion of the sky at one instant and so the 
constraint on a single event with the uncertainty in explosion much greater than the 
orbit of CGRO is weak. 



4. Description of the Concept and Numerical Methods 

We have investigated the question of whether the optical properties of SN 1998bw can be 
explained based on an asymmetric explosion with otherwise normal amounts of kinetic 
energy and total ejecta and 56 Ni mass as opposed to the unprecedentedly large values 
of these quantities required in the "hypernova" models of Iwamoto et al. (1998) and 
Woosely, Eastman, & Schmidt (1998). 

For the initial setup, we use the chemical and density structures of spherical C/O 
cores of Nomoto and Hashimoto (1988). These structures are scaled to adjust the total 
mass of the ejecta. This is an approximation, but the details of the chemical profiles are 
not expected to effect the light curves. 

The explosion models are calculated using a one-dimensional radiation-hydro code 
which includes a detailed nuclear network. The code also simultaneously solves for the 
radiation transport via moment equations. Photon redistribution and thermalization is 
based on detailed NLTE-models. Several hundred frequency groups are used to calculate 
monochromatic light curves, frequency-averaged Eddington factors, and opacity means. 
A Monte Carlo scheme is used for 7-rays. For details, see Hoflich, Wheeler & Thielemann 
(1999) and references therein. 

Aspherical density structures are constructed based on the original spherical density 
distribution. For the simple models presented here, we impose the asymmetry after 
the ejecta has reached the homologous expansion phase. We generate an asymmetric 
configuration by preserving the mass fraction per steradian from the spherical model, 
but imposing a different law of homologous expansion as a function of the angle from 
the equatorial plane. For typical density structures, a higher energy deposition along the 
polar axis results in oblate density structures. Such an energy pattern may be produced 
if jet-like structures are formed during the central core collapse as suggested by Wang 
& Wheeler (1998). In contrast, a prolate density structure would be produced if more 
energy is released in the equatorial region than in the polar direction. For more details, 
see Hoflich, Wang & Wheeler (1999). 

The emissivity along a given isodensity contour in the asymmetric models is taken to 
be constant and equal to the corresponding equivalent velocity-weighted mean layer in 
the spherical model. The bolometric and broad band light curves are constructed by con- 
volving the spherical light curves with the photon redistribution functions, L(O) /L(mean), 
that are calculated by the Monte Carlo code for polarization (Hoflich 1991, Hoflich et al. 
1995). Typical conditions at the photosphere, and therefore the colors, are expected to 
be similar in both the spherical and aspherical models because the energy flux, F(Q), is 
found to be similar both in the spherical and aspherical configuration to within w 40%. 
Stationarity is assumed to calculate the photon redistribution functions since the geom- 
etry does not change during the typical diffusion time scale. We assume implicitly the 
same mean diffusion time scales for both the spherical and aspherical configurations. 
This mostly effects the very early phases of the light curve when the hydrodynamical 
time scales are short. For more details, see Hoflich, Wang & Wheeler (1998). 



5. Results 



We construct models in such a way a way that at day 20 the axis ratio at the photosphere 
is 2. In comparison to the spherical model, the homology expansion parameters are a 
factor of « 2.2 larger along the pole for oblate ellipsoids and a factor of w 1.5 larger for 
prolate ellipsoids along the equator. 

We first calculated aspherical light curves based on the C/O core C021 which gives 
a good representation of the BVRI light curves of the SN Ic 19941 (Iwamoto et al. 
1994). The ejecta mass is 0.8M Q and the explosion energy is Ekin = 10 51 erg. A mass 
of 0.08 Mq of 56 Ni is ejected. This model failed to produce the peak brightness by a 
factor of 2, gave too short a rise time by about 5 days, and shows blue color indices at 
maximum light. 

To boost the total luminosity to the level of observation we increased the amount 
of ejected 56 Ni to 0.2 Mq. This quantity of nickel is still below the estimates for 56 Ni 
of 0.3 M Q in the bright SN II 1992am (Schmidt et al. 1997). As shown in Fig.l, the 
time of maximum light is rather insensitive to asphcricity effects. The need to delay 
the time to maximum and to produce the red color at maximum light suggested the 
need to increase the ejecta mass with an appropriate increase in the kinetic energy to 
provide the observed expansion. We thus computed a series of models with M e j = 2M 
E kin = 2 x 10 51 erg and M Ni = O.2M . 

Asphericity of the amplitude assumed here can change the luminosity over a range 
of roughly 2 magnitudes (Fig. 1). For oblate ellipsoids, the luminosity is enhanced along 
the pole whereas for prolate density structures the enhancement occurs in the equatorial 
direction. Combined with the polarization properties, this provides a clear separation 
between oblate and prolate geometries as P always goes to if the structure is seen 
pole-on and P increases towards lower latitudes (Hoflich 1991). 

Observations of the polarization of SN 1998bw 23 days after the explosion show little 
polarization (< 1%, Patat et al. 1998). By day 58, the intrinsic polarization was reported 
to be 0.5% (Kay et al. 1998). Polarization data on SN Ic is rare, but this value is less 
than seen in some SN Ic and related events (see above). SN 1998bw was also rather 
bright. This combination implies oblate geometries if asymmetry is involved. 

For the comparison between the observed and theoretical light curves as shown in 
Fig. 2, we have used the relative calibration of Woosley, Eastman, & Schmidt (1998) 
for the "bolometric light curve". The broad-band data was obtained from Galama et 
al (1998). To account for the obervations with our asymmtric model, the object must 
be seen from an angle of > 60° from the equator. The same conclusion can be drawn 
independently from the detected, but relatively small linear polarization. This model 
thus suggests that the maximum polarization of SN 1998bw was greater than the 0.5% 
measured by Kay et al. (1998). Overall, the broad-band light curves agree with the data 
within the uncertainties. The intrinsic color excess B-V matches the observations within 
0.1 m and, after the initial rise of ~ 7 days, the agreement in each band is better than 
0.3 m . The main discrepancy with the observations occurs during the initial rise when the 
diffusion time scales are much longer than the expansion time. Under these conditions, 
our approximation for redistribution of the energy of a spherical model breaks down 
since the diffusion time scale is long compared to the hydrodynamical time scale. The 




Fig. 1. Directional dependence of the bolometric light curve for oblate (left) and prolate (right) 
ellipsoids. The luminosity of the corresponding spherical model is shown as L(mean). In ad- 
dition, the instantaneous 7 — ray deposition is shown. P(0°) is the polarization at maximum 
light (P(9) ~ P(0°) x cos 2 9). 



decline after maximum is slightly too steep in the models both in the bolometric and 
broad band light curves. This is likely to be related to the energy generation in the 
envelope by 7-ray deposition or to the change in the escape probability of low energy 
photons. The decline rate immediately after peak can be reduced by increasing the 
amount of radioactive 56 Ni by « 40%. In our models, the escape probability for 7- 
rays increased rapidly between day 20 and 80 from 3% to 50%. An alternative means 
to flatten the light curve is to reduce the increase in escape probability. This can be 
achieved by a modification restricted to the inner layers of the ejecta because the escape 
probability is determined by those layers. Either the expansion velocity of the inner 
layers can be reduced or the density gradient might be steeper. Both are expected for 
strongly aspherical explosions. 

6. Discussion and Conclusions 

The "hypernovae" models for SN 1998bw provide interesting fits to the data, but these 
models have some problems. Although the fit of Iwamoto et al. (1998) of the light curve 



SN98bw: Theory vs. Observation 




time [days] time [days] 

Fig. 2. Comparison of bolometric and broad-band light curves as observed for SN1998bw with 
those of an oblate ellipsoid seen at high angle with respect to the equator assuming a distance 
of 36 Mpc and A v = 0.2 m . 



is excellent with errors < 0.3 m over 40 days, the spectra show absorption lines that are 
too narrow by a factor of 2 to 3 indicating too narrow a range of formation in velocity 
space. This may be related to the high envelope mass. In the lower mass models of 
Woosley, Eastman, & Schmidt (1998), the computed color indices (B-V, V-R, V-I) are 
too red at all epochs. The hypernova models presented to date are spherically symmetric 
and hence make no pretense of accounting for the observed polarization, but this is a 
critical feature and cannot be ignored. 

We have shown that the high apparent luminosity of SN 1998bw may be understood 
within the framework of "classical" SNIc by invoking asymmetry of the ejecta of the 
degree required to account for the polarization. 

Even with the invocation of significant asymmetry, our model for SN 1998bw remains 
at the bright end of the scale for normal SN Ic. We note that the luminosity of SN 1998bw 
may be uncertain by a factor of 2 due to non-Hubble motion within the cluster and 
uncertainties in the Hubble constant and reddening. For a model with an ejected mass 
of 2 Mq, an explosion energy of 2 x 10 51 erg, and a 56 Ni- ejection of 0.2M Q , both the 
bolometric and broad-band light curves are rather well reproduced by an oblate ellipsoid 
with an axis ratio of 2 to 1 which is observed within 30 deg of the symmetry axis. This 



angle for the line of sight is consistent with the low (but still significant) polarization 
observed for SN 1998bw. In a Lagrangian frame, the polar expansion velocity is a factor 
of 2 larger than the mean velocity. This is also in agreement with the rather large 
expansion velocities seen in SN 1998bw. 

Neither this model nor the hypernova models give an obvious explanation of the 7-ray 
burst nor the especially bright radio emission. Woosley, Eastman, & Schmidt (1998) have 
analyzed the possibility of 7-ray bursts in the framework of spherical models. Even with 
their explosion energies of more than 20 foe they showed that the 7-ray burst associated 
with SN 1998bw/GRB 980425 cannot be explained by the acceleration of matter to 
relativistic speeds at shock-breakout. In our picture, the specific energy released in the 
polar region is comparable. We want to stress, however, that the asymmetry in the 
energy distribution in the model presented here is set after homologeous expansion 
has been established. Because the early hydrodynamical evolution will tend to wipe 
out asymmetries, the initial anisotropy in the energy distribution as the explosion is 
initiated after core collapse is expected to be significantly higher. 

We have shown that SN 1998bw may be understood within the framework of "clas- 
sical" core-collapse supernovae rather than by a "hypernova," but the actual model 
parameters must be regarded as uncertain both because of the model assumptions and 
the uncertainty in the observed luminosity. In light of the good fits, however, SN 1998bw 
may indeed be a "hypernova." Continuous measurements of the polarization and the 
velocity of 56 Co lines are critical to unreveal the nature and geometry of this object. If 
the late-time tail tracks the 56 Co decay line, then the ejected 56 Ni mass could be de- 
termined. This might prove the simplest discriminant between the "hypernova" models 
and models based on significant asymmetry. 

There remains the broader issue of the possible connection of supernovae and 7-ray 
bursts. If some of the 7-ray bursts at large redshifts are driven by supernovae, then the 
associated energy flow must be strongly collimated and aspherical core collapse may be 
a natural way to produce such jets. The basic purpose of this work is to underline the 
fact that asymmetries must be taken into account in a complete consideration of core 
collapse supernovae, of SN 1998bw in particular, and, by extension of 7-ray bursts. 
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